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So far the squarks have not been detected at the LHC indicating that they are heavier than a few 
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We study the possibihty that a supersymmetric partner of the top quark, stop, is the next-to-hghtest 
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may be difficult. 
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I. INTRODUCTION 

The next-to-lightest supersymmetric particle (NLSP) 
has a crucial role in the attempts of detecting supersym- 
metry (SUSY) since the lightest supersymmetric particle 
will escape the detectors unnoticed. The Large Hadron 
Collider (LHC) accelerates and collides mainly protons, 
and thus the collision processes are overwhelmed by the 
strong interactions. The superpartner of the top quark, 
stop, can be the lightest colored superpartner, mainly due 
to the splitting of the two stop scalar states amplified by 
the top Yukawa coupling. 

In principle, the masses of the superpartners are free 
parameters in the minimal supersymmetric extension 
of the standard model (MSSM). The large number of 
the parameters parametrize the supersymmetry break- 
ing, which is expected to be spontaneous in a more com- 
plete theory. The explicit SUSY breaking is introduced 
softly so that no quartic divergences re-appear. This re- 
quires inclusion of all the possible breaking terms, which 
are gauge invariant, into the Lagrangian. In models of 
a particular SUSY breaking mechanism, the number of 
parameters may be substantially smaller. Perhaps the 
most studied model is the constrained minimal super- 
symmetric standard model (CMSSM) [1, 2], which we 
also consider here. In CMSSM, supersymmetry is sup- 
posed to be broken spontaneously at the hidden sector 
by the SUSY breaking fields, which do not have gauge 
interactions with the SUSY fields that are, in principle, 
observable to us. 

The breaking fields affect us through the gravitational 
strength interactions, which generate an effective La- 
grangian with SUSY breaking terms. The effective La- 
grangian does not have to be renormalizable anymore. 
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SO it contains terms with couplings to the hidden sector 
fields suppressed by powers of inverse Planck mass. Af- 
ter the hidden sector fields generate vacuum expectation 
values, (renormalizable) soft SUSY breaking terms arise. 
This implies, under certain minimality assumptions, that 
the SUSY scalar terms have common couplings: the com- 
mon scalar mass parameter TOq, a common bilinear pa- 
rameter Bq , as well as a common trilinear parameter Aq , 
at the energy scale where the Lagrangian is established. 
This scale is thought to be the grand unification (GUT) 
scale. 

Since the observations are made at the low energies, 
the renormalization group equations (RGE) must be used 
to calculate the running parameters, like the sparticle 
masses, at the low-energy scale. Even though the scalars 
have equal mass parameters at the GUT scale, this is 
not true at the electroweak (EW) scale anymore. Each 
scalar RGE has terms proportional to the corresponding 
fermionic partner Yukawa coupling. The effect of these 
terms is to decrease the mass parameter in question. Be- 
cause of the large size of the top Yukawa coupling, the 
corresponding mass parameters entering to the stop mass 
matrix tend to decrease most, leading to the situation 
where the stops are the lightest squarks. The terms with 
gauge couplings have an opposite effect, as they appear 
with an opposite sign. Since the slepton mass parameter 
RGEs lack terms with the strong coupling constant, they 
do not increase as much as the squark mass parameters. 
Therefore, both the lighter stop and the lightest slepton, 
stau, are natural candidates for being the lightest super- 
symmetric scalar. 

Supersymmetry by itself does not prevent the introduc- 
tion of baryon and lepton number violating terms, which 
are not present in the standard model (SM) renormaliz- 
able Lagrangian. Such terms have a potential to lead to 
a very fast proton decay. Their presence is prevented by 
so-called i?-parity [3]: SM particles are positively, and 
their superpartners, sparticles, negatively charged under 
this parity. The remarkable consequence is that the light- 



est supersymmetric particle (LSP) must be absolutely 
stable. Therefore, a large number of LSPs may still be 
around from the early universe, which could explain the 
observed dark matter relic density (RD). In particle col- 
liders, sparticles can be produced only in even numbers, 
since the initial state contains only ordinary particles. 
The LSP is thought to be uncharged (in both electric and 
color charges) which means that it may be only weakly 
interacting or that it only has gravitational interactions. 
Therefore, it is expected that it escapes the detectors un- 
noticed; only a missing energy component transverse to 
the coUiding beam, E'^''^^ (MET), is noticed. 

The obvious candidate for the LSP in the MSSM spec- 
trum is the lightest neutralino, Xi, which is the mixture 
of the neutral higgsinos, bino and wino, i.e. the super- 
partners of Higgs, B and M^-bosons, respectively. Other 
candidates are the superpartner of graviton, gravitino, 
and superpartner of a neutrino, sneutrino. In gauge me- 
diated supersymmetry breaking (GMSB) models [4], the 
LSP is always the gravitino, which is thus a natural dark 
matter candidate in these models. In GMSB, the leading 
contribution to the squark masses is proportional to the 
strong coupling constant, being larger than the leading 
contribution to the slepton masses. Thus, it is not ex- 
pected that a strongly interacting particle is the NLSP in 
a minimal GMSB model. The searches at LEP collider 
exclude the possibility of light left-handed sneutrinos as 
the LSP, and masses beyond LEP's reach are ruled out 
by direct detection dark matter searches [5]. Unlike the 
partner of the left-handed neutrino, the partner of the 
right-handed neutrino is a viable dark matter candidate 
[6]. In such an extended model, also signatures of stop 
NLSP have been studied at the colliders [7]. 

In this paper, we study, within GMSSM, the possibil- 
ity that the lighter stop scalar state is the NLSP. This, 
in particular, states that the lightest stop ii is lighter 
than the lighter chargino {xi ), the mixture of charged 
higgsinos and winos^. This also implies that the LSP 
is a neutralino'^. This kind of a scenario requires large 
splitting between the two stop scalars. One consequence 
is that a wide gap opens in the sparticle mass spectrum 
between the scalars. ii is the lonely scalar and close to 
the mass of Xi as a result of the relic density constraint. 
Therefore, stop is supposed to be the sparticle that is pro- 
duced at the LHC in abundance. Stop NLSP in CMSSM 
with nonzero trilinear term has been studied, for exam- 
ple, in Refs. [11-13]. Implications at colliders have been 
studied, e.g., in Refs. [13-20]. In this work, we scan over 
all parameters and thus will comprehensively determine 
the possible regions of stop NLSP. We also map all the re- 
gions where dark matter constraints are met. Stop NLSP 



^ In the anomaly mediated supersymmetry breaking (AMSB) 
models [8] chargino is the NLSP. 

^ The gravitino mass is not fixed in CMSSM, but the gravitino 
LSP with stop NLSP is not allowed [9]. In the so-called NUHM 
models [9, 10] stop can be the NLSP while gravitino is the LSP. 



may also arise in some compressed SUSY models [21, 22] 
or in the context of mirage mediation scenario [23, 24]. 
Light stop scenario is also favored in the view of the b — r 
Yukawa unification [25]. 

In section II, we discuss the conditions for stop being 
the NLSP and the decay modes of the stop NLSP. In 
section III, we list the constraints used in this paper. In 
section IV, we scan over relevant parameters to map the 
possible stop NLSP regions. In section V, we discuss the 
possibilities for detection of stop NLSP at the LHC, and 
we conclude in section VI. 



II. STOP AS AN NLSP 

Stop can be the NLSP when the mixing term of the 
left- and right-labeled scalar states M£^ = v {at sin (^ — 
fiyt cos /3) in the stop mass matrix is large enough. The 
largest mixing occurs when the supersymmetric Higgs 
mass parameter ^ and the trilinear SUSY breaking pa- 
rameter at have opposite signs. Here, yt is the top 
Yukawa coupling, /3 is defined in the relation tan/3 = 
'v^jvx, Vi being the vacuum expectation values of the two 
Higgs doublets obeying the relation v = yjv\ -t- v\. We 
assume here that the trilinear soft parameters are propor- 
tional to the Yukawa couplings, ai — yiAi, so the above 
mixing term can be written as M|^ = mt{At — /iCot/3). 

In addition to the mixing, the renormalization group 
running plays an important role in the determination of 
the NLSP. In the scalar RGEs, the terms proportional to 
the Yukawa couplings (and the scalar masses) decrease 
the soft mass parameters, while the terms proportional 
to the gauge couplings (and gaugino masses) have an 
opposite effect. Because of the largeness of the strong 
coupling constant, the squark mass parameters tend to 
increase more than the slepton mass parameters, even 
though the third-generation squarks have large Yukawa 
couplings. Therefore, a small gaugino mass parameter 
as compared to the scalar mass parameter is preferred in 
order to suppress the strong coupling term in the squark 
RGE, which then may lead to a stop NLSP instead of 
stau NLSP. 

The mass of the NLSP affects the relic density through 
co-annihilations. (The neutralino-stop co-annihilations 
were studied in Ref. [26].) Unless the annihilating LSPs 
are close to a resonance or are light, their annihilation 
usually cannot result in dark matter abundance within 
the observed limits. Rather, the relic density is usually 
too large. Co-annihilations with other particles, how- 
ever, dilute the relic density, and with certain parameters, 
annihilations can even be too effective so that there is 
hardly any dark matter left after the annihilations cease. 
In the case of large (negative) values of At, O(TeV), and 
stop NLSP, the mass difference between stop and the 
neutralino LSP should not be less than 20 GeV or more 
than around 50 GeV in order to obtain the desired relic 
density, as we will show later (in the case of positive /x) . 
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FIG. 1. Stop lifetime [fs] as a function of the NLSP-LSP 
mass difference. The points are the same RD=WMAP points 
as used later in Fig. 9. 



only possible stop decays are 

• 1 1 ^ ux° 

• ii -^ cxI 

• h ^ bff'xl 

Because of the required mass difference, the lifetime of 
stop is short, of the order of 10~^^s (Fig. 1). The 
loop decay ti — >■ cx° is enhanced by a large loga- 
rithm In (AQurp/?7i^), and is likely to be dominant over 
the four-body decay, while smaller Cabibbo-Kobayashi- 
Maskawa (CKM) matrix elements suppress the other loop 
decay [27]. There are also scenarios where the four-body 
decay rate is at the same level or larger than the loop 
decay, for example when squark masses are not unified 
at a large mass scale, resulting only in a small logarith- 
mic enhancement [28]. In section V, we will discuss the 
decay modes of stop when combined with the constraint 
requirements. 



III. CONSTRAINTS 

In order to test the validity of a model, a number of 
constraints must be considered. 



1. The Higgs boson mass and other direct searches 

The nonobservation of Higgs bosons and supersymmet- 
ric partners put stringent constraints on the parameters 
of SUSY models. The 95% C.L. exclusion Hmit for the 
SM Higgs boson mass oimn > 114.4 GeV [29] is applica- 
ble also to the lightest CP-even Higgs boson of the MSSM 
if it is SM-like. This implies that the Higgs coupling to 
the Z boson stays near to the SM value. If the hZZ- 
coupling is reduced, as may happen, e.g., at the large 
tan/3, a general SUSY Higgs mass bound, m^j > 92.8 
GeV, must be applied [29]. 



The light Higgs boson mass in SUSY models is gen- 
erally sensitive to the higher-order corrections. Since 
there is a theoretical uncertainty of some 3 GeV from the 
next-to-next-to-leading order (NNLO) and higher correc- 
tions [30], a limit of rrih > 111 is used in our calcula- 
tions instead of the kinematic bound of 114.4 GeV. The 
SM-like Z-coupling is checked for this bound: A value 
sin^(/3 — aeff) > 0.9 is required for the mh > 111 GeV 
limit in our work, otherwise, a limit mh > 91 GeV is 
used, which also leaves some room for theoretical uncer- 
tainties. Here, cteff is the effective mixing angle between 
neutral Higgs bosons. 

The lightest Higgs boson mass is limited by mz at 
the tree level [31] even for more than two Higgs doublets 
[32]. The Higgs mass depends strongly on the tan/3, and 
a large tan/3 is favored. The radiative corrections are 
crucial for bringing the light Higgs mass up to the ac- 
ceptable level. The largest contributions typically come 
from the top and stop loops. In the decoupling regime 
{ttia >• mz with tan/3 » 1), the leading 1-loop radiative 



[33], where 
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Here, Ms is a common soft SUSY-breaking mass term 
for the third-generation squarks and Xt^b are the mixing 
terms involving ^, tan/3, and the stop/sbottom trilinear 
couplings. 
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(3.2) 



Xt = At- n cot /3, 
Xb = Ah — II tan /3. 



(3.3) 
(3.4) 



Because of the quartic dependence on the top mass, the et 
term is the leading one. In the maximal mixing scenario 
[34], the second and third terms of Eq. (3.1) cancel each 
other, leading to the mixing condition 



Xt - V6Ms 



(3.5) 



In the case of stop NLSP, the Higgs mass constraint seems 
to require nearly maximal mixing, as can be seen later. 

Also, the other collider bounds from LEP, like chargino 
mass limits, and the LHG squark and gluino limits (CMS 
[35], ATLAS [36]) are important constraints. The LEP 
limits have been taken into account, as implemented in 
the program micrOmcgas (v.2.4.R) [37]. The new LHG 
limits were also considered, where applicable. 



2. The B -^ TV decay 

The final states of the decays of the type B ^ Iv con- 
sist purely on leptonic states, and the hadronic uncer- 
tainties are present only in the B meson decay constant 



Jb- Only the r channel has been observed so far. The 
SM expected value for the branching ratio, when using 
the value of \Vub\ given by the UTFit collaboration [38], 
is predicted to be [39] 



BR(B+ -^ T+iyr)sM = (0.80 ± 0.15) x 10" 



(3.6) 



Recently, the experimental measurements from the B- 
factories have improved significantly, the most recent 
world average measurement being [40] 



BR{B+ -^ r+i/^)exp = (1.68 ± 0.31) x 10" 



(3.7) 



The new physics (NP) contribution to the branching ratio 
can be quantified by defining a ratio [41] 



R 



NP 



BR{B+ -^ T+Ur) 



SM+NP 



BR(B+ -^T+Vr) 



SM 



(3.8) 



where the subscript SM+NP denotes the branching ratio 
of the NP scenario, including the SM contribution. The 
95% C.L. allowed range for fl^P is then [39] 



0.99 < i?^P < 3.19. 



(3.9) 



With higher-order corrections, the formula for the ratio 

is [41, 42] 



i??J=(l- 



tan^ 



1 + eo tan/3 M^± 



(3.10) 



where eo contains all the higher-order corrections and 
ruB is the B-meson mass. The branching ratio depends 
strongly on the charged Higgs mass and tan/3, and typi- 
cally the constraint (3.9) prefers small values of tan/3 in 
order not to decrease the ratio too much below the lower 
limit (unless M^± is very light). A large charged Higgs 
mass decreases the new physics contributions in general. 



This is a 3.4 a deviation. The SM prediction is largely 
dominated by leptonic QED processes, though careful 
calculation of hadronic and electroweak contributions is 
also necessary due to high precision of the experimental 
measurement. The hadronic processes, vacuum polariza- 
tion above all, have more than an order of magnitude 
larger contribution to the magnetic moment compared 
to the electroweak processes. 

Purely supersymmetric contribution of MSSM to a^ 
is proportional to tan/3sign(/i)/Af|ugY, and for large 
enough tan/3 and not too heavy supersymmetric parti- 
cles, it can be larger than the electroweak contribution. 
For positive values of /x, the MSSM can provide the so- 
lution to the discrepancy. For negative values of /i, the 
new physics contributions drive the gap even wider. Be- 
cause of this, a negative sign for ji in the MSSM is usually 
considered to be disfavored [46]. 

As pointed out in [39], this constraint is complemen- 
tary to the B ^ TV decay constraint, and taken together, 
they rule out large areas of the parameter space. 

The theoretical calculation of a^ is known to be diffi- 
cult because of the hadronic contributions and nonper- 
turbative effects involved, see e.g. [47-49] and references 
therein. Although there has been impressive improve- 
ments in the calculation, frequently the a^ constraint is 
not used in determining the excluded parameter space. 
We will take a similar attitude here, but comment on 
flp when appropriate. When referring to this constraint, 
the following 95% C.L. limits, which include theoretical 
uncertainties, are used [50]: 



1.15 X lO"-* < Aa 



SM+NP 



< 4.75 X 10" 



(3.14) 



Our acceptable parameter points have Aaf^^'^^^^ below 
this range. 



4- The 6 — >■ S7 branching ratio 



3. The {g — 2) o/ muon 

The anomalous magnetic moment of muon, a^ = {g — 
2)/2 has been measured quite precisely a decade ago. 
The measured value [43] for a^ is 



a^fP = (11659208.0 ±6.3) x 10"^ 
The recent SM prediction for this is [44] 



(3.11) 



a™ = (11659178.5 ± 6.1) x 10"^°. (3.12) 

leading to a discrepancy between the SM and experi- 
ment'^. 



Aa^ = (29.5 ±8.8) x 10 



-10 



(3.13) 



The present experimental value by the Heavy Flavor 
Averaging Group (HFAG) is [51] 

BR{B -+ Xs-f) = (355 ± 24 ± 9) X 10"^. 

In our constraint, the theoretical uncertainties are in- 
cluded as well: in the SM at the NNLO QCD level the 
uncertainty can be estimated to be 23 x 10^^ [52], in the 
MSSM, the theoretical uncertainty is estimated to be ad- 
ditional 5% (we take 15 x 10^^) [53]. Combining all these 
gives (at 2 a) 

BR{B -^ X,-i) = (355 ± 142) x 10"^ (3.15) 

The b -^ s"f constraint is sensitive to the sign of fi [54], 
preferring the positive value. 

5. Other constraints 



^ In a very recent paper [45] it was argued that the experimental 
and SM value actually agree. 



The Bs — ^ M"*"/^ branching ratio can also be used as a 
constraint for new physics. The most recent experimental 



upper limit by the CDF collaboration is [55] 

BR{Bs -^ M+/"") < 4.3 X 10"* (95% C.L.) (3.16) 

Including the theoretical uncertainty from fss = 238.8 ± 
9.5 MeV [56], we find a conservative upper limit of 



BR{Bs -^ fi'^fi^) < 5.0 X 10" 



(3.17) 



In practice, this constraint is usually not the limiting fac- 
tor, and in our calculations it is never the main constraint 
in otherwise good points. Therefore, we do not comment 
on this constraint further. 



the 4-dimensional parameter volume of the three GUT 
scale parameters and tan/3 and calculated the low en- 
ergy particle spectrum and the constraints for each point. 
The sign of /i was chosen to be positive, partly guided by 
the preferences for some of the constraints of Sec. III. It 
should be noted that we restricted ourselves to the se- 
lection of real-valued parameters, which implies that no 
additional CP-violation is introduced by the soft SUSY 
breaking terms. 

The particle spectrum was calculated using SOFT- 
SUSY (v. 3. 1.7) [60], and the relic density and constraints 
using micrOmegas (v.2.4.R) [37]. Top pole mass of 
nit — 173.3 GeV was used throughout this study. 



6. Relic density 

Stable, (color and charge) neutral LSP provides a con- 
venient candidate to explain the dark matter abundance 
observed by the Wilkinson Microwave Anisotropy Probe 
(WMAP) satellite mission. With the combined data from 
7-year WMAP results, BAO (Baryon Acoustic Oscilla- 
tions) in the distribution of galaxies, and observation of 
Hubble constant, the density of cold dark matter in the 
universe is determined to be [57] flch^ — 0.1126±0.0036. 
If 10 % theoretical uncertainty is added [58], we find the 
preferred WMAP range of 



0.0941 < n^h^ < 0.1311 



(3.18) 



at 2 a level. CMSSM can provide a dark matter density 
that is within the WMAP limits, at least in some parts 
of the parameter space. However, dark matter may be 
of nonsupersymmetric origin, so here we have used only 
the upper bound as a real constraint, unless otherwise 
indicated. 



IV. STOP NLSP PARAMETER SPACE 

A. Method 

The CMSSM parameter space is in practice four di- 
mensional (?7io, A'/i/2, ^Oi tan/3), with the sign ambigu- 
ity of the /^-parameter doubling it. The actual value 
of the /i-parameter is calculated from the radiative elec- 
troweak symmetry breaking (rEWSB) constraint condi- 
tion [31, 59], which leaves the sign of /^ as a free pa- 
rameter. Ml/2 is the common gaugino mass, defined 
at the GUT scale by the boundary condition A/1/2 — 
Ml = M2 = M3, where Mi,2,3 are the [/(I), SU{2) and 
SU{3) gaugino masses, respectively. Limiting oneself to 
2-dimensional parameter scans, on the one hand, reduces 
the effect of the other parameters to discrete examples, 
and on the other hand, creates easily an illusion of very 
restricted parameter space volume. Therefore, it is useful 
to actually explore the full parameter space. In order to 
find the parameter points, where stop is the NLSP and 
the experimental constraints are not violated, we scanned 



B. Stop NLSP regions 

Figure 2 shows a 2-dimensional projection of the 
scanned 4-dimensional parameter space, where the 
ranges for parameters are (in [GeV], fi > 0): 

50 < Ml/2 < 1490|20; 

50 < mo < 1990|20; 

-5000< Aq <5000|100, 

tan^e {2...7|1;10...50|10}. 

The number after "|" denotes the stepping. Each point 
shows the identity of one of the possibly many potential 
NLSPs which results from varying the remaining two free 
parameters (tan/3 and toq in this case). The representa- 
tive NLSP identity is chosen by keeping in mind that we 
would like to find all the points where stop can be the 
NLSP. Therefore, if a stop NLSP is found, it is chosen. 
Otherwise a stau, xt or X2 NLSP is chosen, in that or- 
der of preference. Thus, an area labeled with fi does not 
allow stop NLSPs for the scanned parameters, but may 
contain, e.g., X2 NLSP. The rge-denoted (yellow) area 
shows the space where no combination of the parame- 
ters provides a good solution to the spectrum calculation 
(tachyons, no rEWSB . . . ). A large area of stop NLSP 
is found when the Ao-parameter is nonzero. With our 
choice of positive fj,, the negative Ao-parameter leads to 
larger mixing in stop mass matrix, see Eq. (3.3). There- 
fore, for a symmetric range, a larger number of potential 
stop NLSP points is found with negative Ag-parameters. 
A large cot /3 is favored because of the enhancing effect 
in the stop mixing on the one hand, and because of the 
suppressing effect in the stau mixing on the other hand. 
Since tan/3 is paired with the /i-parameter, which is de- 
termined from the rEWSB condition, general upper or 
lower limits for tan/3 are very involved. From our numer- 
ical calculations, however, we learn that the stop NLSP 
exists up to tan/3 — 50, but for tan/3 = 55, stau is the 
NLSP. For a large value of M1/2 , the upper value for tan /3 
drops to about 35 (stau being the NLSP otherwise). To 
have a better view of the stop NLSP points, we look at 
the negative ^o'P^-ra-nieters in more detail. We anticipate 
this also to alleviate the lightest Higgs boson mass limit 
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FIG. 2. The potential stop NLSP parameter space in dark 
(brown) on (Mi/2,^0) plane. No constraints are applied yet. 



constraint, as discussed in Sec. Ill 1. Some numerically 
unstable points may give erroneously large vrih, hence, we 
set also an upper limit for rrih to be 140 GeV [33, 61, 62], 
which cuts the unreliable points from the calculations. 

In Fig. 3, the negative Ap-parameter space is plotted 
with respect to M1/2 as a 2-diniensional projection of the 
4-dimensional parameter space showing the NLSP map. 
Parameters for subsequent figures are (in [GeV], /i > 0): 

100 < A/i/2 < 1800|20; 

1000 < mo < 2848|66; 

- 8040 < Ao < -2200|80; 

tan/3e {2. 5.. .20(0. 5; 25. ..55(5} 

A large space with stop NLSP is found. However, many 
of the potential stop-NLSP points conflict with the con- 
straints specified earlier. The requirement of the neu- 
tralino relic density to fully explain the observed WMAP 
cold dark matter density is quite a restrictive constraint. 
This constraint may be weakened to the requirement that 
neutralino dark matter relic density should not exceed 
the upper WMAP limit for relic density. This, however, 
means that we cannot explain the dark matter puzzle 
with the SUSY model. Figure 3 shows, on top of the 
NLSP map, also scattered (red) points, where stop is 
the NLSP and the neutralino relic density is below the 
WMAP upper limit. The white subset of these points 
shows the parameters for which, in addition, the other 
constraints of Sec. Ill are not violated (however, exclud- 
ing the muon magnetic moment, which is below the range 
(3.14)). In other words, these points are not excluded 
by the constraints discussed and exhibit the stop NLSP. 
When examined one by one, the BR(i3+ — > t^Vt) con- 
straint limits the points only near the rge-denoted area, 
which is true also for the 6 — > 57 constraint. In contrast, 
the possible discrepancy in {g — 2)^ cannot be explained 
with these models, since the calculated SUSY contribu- 
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FIG. 3. The potential stop NLSP area in (Mi/2,^0) plane. 
Zoomed from Fig. 2 to an interesting area. The red points 
are the stop NLSP points with relic density below the WMAP 
upper limit (which we denote from now on as RD < WMAP) , 
and the white subset of the points also obey the constraints 
of Sec. III. See the text. 



tion here is always below the range of Eq. (3.14). When 
applying all the constraints (except {g — 2)^) simultane- 
ously, the white point area appears. It should be noted 
that there exists a subset of the above points, which 
nicely match the WMAP constraint, thus explaining also 
the amount of dark matter in the universe. Moreover, it 
matches the shown white dot area, just being sparser due 
to the limited number of scanned points. Therefore, it is 
possible to find points where stop is the NLSP, neutralino 
relic density matches the WMAP observation, and the 
collider constraints are fulfilled. This leaves still the pos- 
sible deviation of the observed {g—2)^ unexplained. Even 
though the [g — 2)^ value could be increased by consid- 
ering the points with large tan/3, this conflicts with the 
BR(i3+ — >• T^Vr) constraint. In addition, \Aq\ should 
not be larger than 2-3 TeV for the [g — 2)^ constraint to 
be fulfilled, [g — 2)^ disfavors the stop-NLSP scenario. 
One may wonder if the large \Aq\ values depicted can be 
physical, since typically charge and color-breaking min- 
ima occur for \Aq\ much larger than rtiQ [63]. However, 
if the Universe is in a false vacuum, the tunneling prob- 
ability to the real minimum can be very small [64]. We 
have checked for some examples with |^o| ~ ^ttiq that 
even if CCB minimum exists, the tunneling time scale is 
longer than the age of the Universe. 

Figure 4 shows the equivalent of Fig. 3, but now the 
points are projected to (mo, Ao)-plane. Figure 5 shows 
the projection to the (mo, Mi/2)-plane. In Video 1, the 
tan/3 = 10 slice of Fig. 5 is plotted. The points of Higgs 
maximal mixing (dark green) are distinguished from the 
other points that obey the relic density upper-limit con- 
straint (red dots) and the points fulfilling also the other 
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constraints (white dots). The green points he within 5% 
range from the optimal mixing value of Eq. (3.5). (See 
the animation of the evolution of the trilinear coupling 
from [65].) The maximal mixing coincides well with the 
good points, suggesting that the Higgs mass limit is, in 
fact, quite a severe constraint. The large area of suffi- 
ciently low relic density points (red) is generated, when 
the RGE conditions push the Afi/2 parameter larger as 
the absolute value of the trilinear-parameter increases. 
This can be clearly seen in the video. 

Figures 3-5 are different projections of the same pa- 
rameter volume. It is interesting to see how the viewpoint 
exaggerates or understates certain aspects of the param- 
eter space. For example, the area with no good spectrum 
(denoted by rge) seems to be quite large in Fig. 3, whereas 
in Fig. 5 it is hardly visible. An important message to be 
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Video 1. Points of optimal mixing (green) for stop NLSP 
points, tan/3 = 10, RD < WMAP. The scan is somewhat 
denser in this plot than in Fig. 5. 



learned from this is that 2-dimensional extracts from a 
multidimensional parameter space may give a false feel- 
ing of the parameter space being very restricted. 



C. Masses 

Searches for squarks and gluinos at the LHC disfavor 
a low mass region [35, 36]. If only stop is light, the an- 
ticipated exclusion limit would be lower than in the case 
of degenerate squarks and gluino due to smaller produc- 
tion cross section, and partly also due to tt background. 
Also, proximity of the stop NLSP to the LSP might result 
in the supersymmetric events failing the missing energy 
cuts. A CMS search for jets-|-MET resulted in Nmax = 13 
events at 95% confidence level for an integrated luminos- 
ity £ = 35 pb~^. The upper bound is related to the total 
SUSY production cross section by iVmax = C'CcTmax- If we 
assume the total efficiency e to be 25% {e.g. see the dis- 
cussion on efficiency in [35]), we can estimate the upper 
bound on the cross section to be cTma 



= 1.5 pb. Esti- 
mating that about one picobarn production cross section 
is needed to exclude a certain sparticle, we get for the 
lonely stop NLSP an anticipated lower mass limit of 300 
GeV from tt* production using Prospino [66]. However, 
stop NLSP may not be easy to detect at the LHC, as will 
be discussed in the next section. 

Interestingly, a small mass difference between stop and 
the LSP, and also a wide gap between the lighter stop and 
the other squarks are exactly what we find in the case 
of stop NLSP. Large splitting in the stop sector pushes 
lighter stop down to be a lonely SUSY scalar. The other 
scalars, except the light Higgs, are much heavier, with 
masses typically above 1 TeV. This is because the stop 
NLSP prefers large tuq as compared to M1/2 (about twice 
the Ml/2 value in the following). In Fig. 6, a typical mass 



> 

o 
O 



"touu 


A 


' 


' 


jKh^ 


4000 


9 \ . 


._^ 


>^.... 


■■■-■-^■^i' 


3500 










^l \ K AH-'- 


y'("" 






3000 


Px -, ;/ ,.,-v.' 


[;...-■ 




:-- 




fi rvP'^'^--'-'' 


.-■\ \ 


, .'. ^■^ 


../*■■■■" "■^" 


2500 


>il 1X2 ...^^...^..^.....^l ..p.. 


: /■ 






2000 


^tJJ/'''..M^^''' '' 








1500 


\-j|fv/ *''■ tS^'-' 






^;j_: 




-%',^-^\ 




._..-— 


■■"■ 


1000 


A^f.»^*> i ,^,- 


-^■':''' 




y ■ 


500 


-■-iIS$----^^ 


i 




i ■ 



200 



300 



400 



500 



600 



700 



800 



FIG. 6. Mass spectrum of a set of good stop NLSP points as a 
function of the stop mass. In the key, particles are hsted in de- 
scending mass order found at mj^ = 800 GeV (RD=WMAP). 
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RD < WMAP. Vertical line reminds about the assumed LHC 
stop mass limit. 
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spectrum is plotted as a function of m^^ . Only one stop 
and x? are light, with masses close to each other. Each 
point in the plot is selected from the white points of the 
scan in Sec. IV B. The stop mass follows quite faithfully 
the Ml/2 parameter, but the other parameters do not in- 
crease monotonically with the stop mass for the selected 
points. Therefore the neighboring points are not neces- 
sarily next to each other in the 4-dimensional parameter 
space. They, however, obey all the constraints except 
a^i, and were also required to explain the amount of dark 
matter (relic density is within the limits of Eq. (3.18), 
which we denote as RD=WMAP). In Fig. 7, we show the 
ratios of several parameters for the points in Fig. 6 (and 
consequently, for Figs. 10-13). In particular, it turns out 
that for these points, |^ol < 3.5too always. The favored 
value of tan/3 ^ 15. It is also clear that for heavier stops. 



Too and M1/2 come closer to each other, while \Aq\ value 
increases compared to toq. 

Figure 8 shows the allowed stop NLSP points in the 
(to.( ,TO.g) plane for RD < WMAP. These points are the 
white points of the scan of Sec. IV B. There is a strong 
correlation with stop and gluino masses. If the stop 
NLSP mass is required to exceed 300 GeV, it would 
imply that the gluino mass is rather heavy. This is 
easily understood, since the LSP is a bino in a large 
part of the CMSSM parameter space. The relic den- 
sity constraint requires the stop mass therefore to be 
rather close to the x? mass, so that co-annihilations 
are able to dilute the excess LSP density. The gaug- 
ino mass relation Mi : M2 : M^ ~ 1 : 1.9 : 6.2 at 
the EW-scale [24] then predicts the gluino mass to be 
about 5-6 times the stop mass {nig sa 6.2 x Mi, and 
TOj^ = rriLSP + 10% = 1.10 X Ml, hence mg/mi_^ w 5.5). 

In Fig. 9, the mass difference of NLSP-stop and Xi is 
plotted with respect to the stop mass for the same al- 
lowed stop NLSP points. The mass difference seems to 
be typically below approximately 50 GeV. This is a con- 
sequence of the relic density constraint, which requires 
effective co-annihilations for the bino LSP. (The exact 
WMAP preferred region saturates the upper edge of the 
mass range, see Fig. 9.) Another implication is that the 
stop decay channels are limited to the ones explained in 
Sec. II. 



V. STOP NLSP AT THE LHC 

In the stop NLSP case, it is quite possible that most 
of the supersymmetric partners are relatively heavy as 
discussed in the previous section. Discovering stop has 
been studied in detail in several works, e.g., in [13-20, 
67, 68]. Here, we will shortly discuss the cross sections 
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FIG. 9. Mass difference of the two lightest SUSY particles 
wrt. ruf^, for black dots RD < WMAP and for light blue 
dots RD=WMAP. Vertical line reminds about the assumed 
LHC stop mass limit. 
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TABLE I. Cross sections [fb] at the LHC for M1/2 = 620 
GeV, mo = 1528 GeV, Ao = -3880 GeV, tan/3 = 14, and 
sign(/i) = +1, corresponding to raj: = 304 GeV of Fig. 6 
(Prospino2, NLO). 



and decay modes of stop NLSP at the LHC, when the 
constraints discussed in previous sections are fulfilled. 

Typically, the tit^ cross section is very large compared 
to the other production mechanisms via cascade decays 
of other sparticles. For example, in Table I, the next- 
to-leading order LHC cross sections for several squark 
and gaugino production channels are calculated for one 
acceptable parameter point corresponding to mi_^ — 304 
GeV of Fig. 6. 

In Fig. 10, we have plotted the dominant decay modes 
of the lightest stop as a function of its mass for the 
points in Fig. 6. The dominant decay mode is through 
x\c channel, but the importance of bff'xl channel in- 
creases with increasing stop mass and the mass differ- 
ence to the LSP (Fig. 9). In Fig. 11, the branching ratios 
of gluino are plotted as a function of its mass for the 
same points. Gluino decays dominantly to a stop-top 
pair. The charge conjugated mode is also included in 
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FIG. 10. Branching ratios of lightest stop for the collection 
of allowed points, RD = WMAP. 
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FIG. 11. Branching ratios of gluino for the collection of al- 
lowed points, RD = WMAP. (*) indicates that the charge 
conjugate state is also included. 



the plotted value. For the evaluation of branching ratios 
and the decays of the supersymmetric particles, we have 
used SUSY-HIT (v. 1.3 with SDECAY vl.3b/HDECAY 
v3.4) [69]. It should be noted, though, that the ap- 
proximate result for the ti — > Xi^ decay used in [69] 
has to be taken with care in case the minimal flavor 
violation scale is not large and has to be reanalyzed 
more carefully [70]. The cross sections have been cal- 
culated to next-to-leading order (NLO) using Prospino2 
(v."on_the_web_ll_3_10")[66]. 

In Fig. 12, the LHC 14 TeV cross section for the like- 
sign top process pp ^ gg ^ XiXiCctt/xiXi^cc [pb] has 
been plotted as a function of stop mass up to 400 GeV, 
where the cross section is already hopelessly small. The 
points correspond to the previously used stop NLSP set. 
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FIG. 12. Cross section of the process pp ^t gg ^r XiXiCctt(*) 
[pb], for collection of allowed points, RD = WMAP. The final 
state includes both combinations of like-sign tops. 



which passes the constraints (points in Figs. 6, 10, 11). 
The plotted cross section contains contributions from 
both charge conjugation final states. For mi = 300 GeV, 
the cross section is around 10 fb but decreases fast with 
increasing stop mass, which is caused by the declining 
gluino production cross section. The background of like- 
sign decays can be removed [18], which makes the process 
interesting for probing light stops at the LHC. 

The stop pair production cross section is clearly the 
largest production channel. However, for a light stop, 
the loop decay dominates, and the final state of almost 
back-to-back neutralinos and two soft c-jets is not exper- 
imentally promising. For heavier stops, the four-body 
branching ratio competes with the loop decay one. In 
Fig. 13, the LHC 14 TeV cross section for the process 
pp — ^ 1x1*1 — 7> (xi6//')(x5c)-fcharge conjugated final 
state [pb] has been plotted for the same stop NLSP set, 
which passes the constraints (the middle curve). The 
cross section remains reasonably large, 0(10 fb), even 
for mi ~ 800 GeV. When the fermions //' are iv from 
W, the signal would be a 6-jet and charged lepton in one 
hemisphere and missing energy from neutralinos, neutri- 
nos, and a soft c-jet in the other hemisphere. A detailed 
signal analysis is not the purpose of this work, but this 
signature may be possible. The cross section where both 
stops decay to four particles is for all studied stop masses 
around 0(10 fb), although the variation between neigh- 
boring points may be large. 



VI. SUMMARY AND DISCUSSION 

We have studied the interesting possibility of stop be- 
ing the next-to-lightest supersymmetric particle within 
the CMSSM scenario. Large mixing in the stop mass 
matrix can result in a (relatively) light stop squark, and 
in some cases, it can be the NLSP. Typically this prefers 
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FIG. 13. Cross section of the process pp — > titi — >■ XiXiCC 
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the middle, including the charge conjugated final state), and 
the process pp — >• iiil — i> {xibfif")ixibfj.f') [pb] (the lowest), 
for collection of allowed points, RD = WMAP. 



a large nonzero value for the trilinear mass parameter 
At- In addition, a large /^-parameter naturally splits 
the stop pair, and may also result in a light mi , unless 
there is cancellation between fi and At . Large parameter 
spaces with the lighter stop as an NLSP can be found, 
which simultaneously fulfill various experimental con- 
straints. Within these regions, the neutralino LSP may 
explain the dark matter problem. The BR(i3+ — > t^i/t), 
BR{B — >■ Xsj) and BR{Bs -^ fi^fi") constraints agree 
with the experimental limits, but we have not found 
points where the possibly large observed discrepancy of 
{g — 2)^( from the SM expected value could be explained. 
The CP-even Higgs boson lower mass limit is a severe 
constraint, which many times rules otherwise good pa- 
rameter spaces out. 

The LHC is currently pushing squark and gluino lower 
mass limits to higher values. The nonobservation of spar- 
ticles may fit to the scenario, where the NLSP is a rather 
heavy squark. In CMSSM, this would in most cases be 
the lighter stop squark. We have argued that observ- 
ing the top signal resulting from this scenario is quite 
challenging. The like-sign top channel, that relies on the 
Majorana character of the gluinos, requires reasonably 
large production cross section of gluinos at the LHC. In 
the stop NLSP scenario with, let us say, heavier than 
350 GeV stop, the gluino is very heavy, about five times 
that, by the WMAP constraint and gaugino mass ratios. 
Therefore, the gluino production will be scant, resulting 
in the like-sign top signal to lie in a subfemtobarn range. 
Because of the relative lightness of stop with respect to 
the other SUSY scalars, their direct pair production dom- 
inates. In the case of heavy stop NLSP, the production 
and decay chain pp — >■ iii^ — ^ biv + I^t may be possi- 
ble (see also [68]). Alternatively to the stop production 
alone, signals where other particles are used in order to 
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discover stops can be utilized, as studied, e.g., in [15] for 
light stops. 
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